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We propose an experimental protocol to directly observe the Kondo effect by scattering ultracold 
atoms. We propose using an optical Feshbach resonance to engineer Kondo-type spin-dependent 
interactions in a system with ultracold ®Li and ®^Rb gases. We calculate the momentum transferred 
from the ®^Rb gas to the ®Li gas in a scattering experiment and show that it has a logarithmically 
enhanced temperature dependence, characteristic of the Kondo effect and analogous to the resistivity 
of alloys with magnetic impurities. Experimentally detecting this enhancement will give a different 
perspective on the Kondo effect, and allow us to explore a rich variety of problems such as the 
Kondo lattice problem and heavy-fermion systems. 
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I. INTRODUCTION 

Ultracold atomic gases provide a platform to engi¬ 
neer model Hamiltonians relevant for condensed matter 
physics phenomena. One such intriguing phenomenon 
is the Kondo effect [Ei- In this paper we propose an 
experimental protocol to engineer and measure the scat¬ 
tering properties of Kondo-like interactions between ul¬ 
tracold atoms. Such an experiment would give a new 
perspective on an iconic problem. 

The Kondo effect is a transport anomaly that arises 
when itinerant electrons have spin-dependent interac¬ 
tions with magnetic impurities. The source of the phe¬ 
nomenon is a spin-singlet many-body bound state formed 
between the Fermi sea and an impurity. This bound state 
leads to resonant scattering of itinerant electrons off the 
screened impurities. As the temperature is lowered, this 
resonant scattering dominates over other scattering pro¬ 
cesses and leads to a characteristic logarithmic tempera¬ 
ture dependence of the resistivity of the material. When 
the interactions between the electrons and the impurity 
are spin independent, no such bound state is formed, and 
the scattering is not enhanced. 

Despite intense research, some questions about the 
Kondo effect remain unresolved and some of the key the¬ 
oretical predictions have never been directly seen. For 
example, the electron cloud which screens the spin on the 
impurity has never directly been imaged il. More im¬ 
portantly the analogous problem with an array of inter¬ 
acting impurities (the Kondo lattice) has aspects which 
are not well understood Q- Exploring the Kondo lat¬ 
tice problem is of paramount importance to the under- 
standiM of heavy fermion systems and quantum critical¬ 
ity 0 ■ 

In this paper we propose using cold atoms to directly 
observe enhanced Kondo scattering. We envision a sys¬ 
tem consisting of a spin-1/2 Fermi gas and a dilute Bose 
gas with spin S, where bosonic atoms play the role of 
magnetic impurities and fermionic atoms play the role of 
electrons. To strengthen the analogy with immobile spin 
impurities in the Kondo model, we consider bosons which 


are much heavier than the fermions. Fermion-boson pairs 
such as ®Li-®^Rb ^Li-®®Rb or ®Li-^^^Cs are good candi¬ 
dates with large mass ratios. Alkaline-earth-metal and 
rare-earth atoms are also promising. 

We consider a rotationally symmetric interaction be¬ 
tween the ultracold atoms, which includes both density- 
density and spin-dependent interactions. We present an 
experimental protocol to produce such an interaction us¬ 
ing an optical Feshbach resonance. For this general in¬ 
teraction, we calculate that the scattering cross section 
is strongly enhanced by the Kondo effect. We propose 
directly measuring this enhancement by launching the 
Bose gas into the Fermi gas with a small velocity. One 
would then measure the momentum transferred to the 
Fermi gas. A number of related experiments have been 
used to probe atomic scattering in the past pTl - [T3l| . We 
show that at temperatures smaller than the Fermi tem¬ 
perature, the final momentum of the Fermi gas varies 
logarithmically with temperature, analogous to the resis¬ 
tance of electrons in an alloy with magnetic impurities. 
The temperature dependence of the transferred momen¬ 
tum, depicted in Fig. [1] has a minimum which is a sig¬ 
nature of the Kondo effect, and this minimum can be 
detected at experimentally accessible temperatures. Al¬ 
ternatively, the enhanced scattering could be seen in the 
dam ping of collective modes of the atomic clouds in a 
trap [l^ . 

This paper is organized as follows. In Sec. [TTlwe in¬ 
troduce our atomic system and the model we consider. 
In Sec. m we explain how an optical Feshbach reso¬ 
nance can be used to produce the interactions consid¬ 
ered in our model. In Sec. m we calculate the mo¬ 
mentum exchanged in a scattering experiment between 
atomic clouds. We calculate the momentum transferred 
as a function of temperature perturbatively up to third 
order in the interaction strength. We explicitly describe 
all parts of our calculation in the appendix. We summa¬ 
rize in Sec. El 
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FIG. 1: Temperature dependence of the momentum P trans¬ 
ferred from bosons to fermions in a scattering experiment with 
photoinduced interactions. Both P and temperature have 
been rescaled to dimensionless quantities. Po denotes the 
momentum transferred to the Fermi gas at zero temperature 
when the interactions are spin independent {ga = 0). (Solid 
line) Spin-dependent interactions between spin-1/2 fermions 
and spin-1 bosons with Qg = ^Qn = ’ (Dashed line) 

Spin-independent interactions (gg = 0 and = —3 x ). 

The minimum in P is a signature of the Kondo effect, and may 
be detected experimentally. In Sec. IIIII we estimate experi¬ 
mental parameters to achieve the interaction strength used 
here. Inset shows a cartoon of the collision. 


II. MODEL 


In this section we describe our model. In Sec. m we 
describe how to experimentally implement our model. 

We build our system out of spin-1/2 fermions and 
spin-S' bosons. In our implementation these will be hy- 

perfine spins. We let the operators dl^ and bl^ create 
fermionic and bosonic atoms at position r and spin pro¬ 
jection a =t, 4. or /X = —S', ..,S along the z-axis. Their 
Fourier transforms. 






( 1 ) 


create particles in momentum eigenstates. Above, V is 
the volume of the system. 

We explore a model with a Hamiltonian H = Ffo+^int- 
The first term models the kinetic energy of the fermions 
and bosons, 


For the interactions modeled by ISmt, we consider a 
generic form of local spherically symmetric pairwise 
Bose-Fermi interactions. Since the fermions have spin- 
1/2, the most general such interaction has the form, 

Hint = I d r ^ ^ d\.Q^dr^yi.^ry “b ■ 

0.(3 

. 

We denote the vector of spin matrices for the fermions 
and bosons by and and 5 refers to the Kro- 
necker delta function. It is important to note that Hiut 
contains terms where a ^ jd and ^ ^ v. This encodes 
the fact that the atoms exchange spin when they collide. 
We point out that spherical symmetry of the Hamilto¬ 
nian is not a necessary feature to observe Kondo physics. 
Any Hamiltonian which allows spin exchange processes 
at third order of interaction strength would produce an 
enhanced scattering cross section at low temperatures. 
We restrict ourselves to interactions modeled by Eq. ©, 
and we show in Sec. imi that this has a simple experi¬ 
mental realization. 

It is useful to rewrite Hint in momentum space as 

^int = J d^k J d^p J d^q Y, (4) 

^k+q,a^k+p,/3bl.-g^fjbk-p,v ‘ . 

Our model in I^. (|T|) differs from the one in the spin-iS 
Kondo model [i| in two respects. The bosonic atoms, 
which play the role of impurities, are mobile. Due to 
their large mass however, the recoil of the bosonic atoms 
can be neglected, and formally the physics is equivalent 
to that of immobile spin impurities. In addition to the 
regular spin-S' Kondo-like interaction, Eq. contains 
a density-density interaction. We show that in spite 
of such an additional interaction term, the momentum 
transferred to the Fermi gas in a scattering experiment 
still has a minimum at a certain temperature, albeit at a 
lower temperature than the case with no density-density 
interaction. 

The interaction we have considered in Eq. m does 
not occur in typical cold atom experiments in which in¬ 
teraction strengths are tuned using a magnetic Feshbach 
resonance. In a typical magnetic Feshbach resonance, 
spin-exchange collisions are off-resonance and will not be 
observed. In the following section we propose using an 
optical Feshbach resonance to produce the interaction in 
Eq. dSD . 

III. AN EXPERIMENTAL SETUP 


Ho = 



^(Cfc - p)dl^dka + Y ^kbY^kp. 


_ ^ 

“ 2ma ’ “ 2Mb ' 


( 2 ) 


In this section we describe our proposal to experimen¬ 
tally implement the model introduced in Sec. El using 
®Li and ®^Rb atoms as our itinerant fermions and spin 
impurities. As we will show, producing a strong inter¬ 
action between ®Li and ®^Rb using an optical Feshbach 
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resonance requires a large matrix element for photoasso- 
ciation. Experiments [l^ show ^Li and ®®Rb to have the 
highest photoassociation rate coefficient among all the 
bialkali metal combinations. We expect their isotopes 
®Li and ®^Rb to have similar photoassociation rates, and 
we chose ®Li and ®^Rb in our proposal to produce the 
Kondo model because they are readily available in ultra¬ 
cold atomic experiments. The ®Li and ®^Rb atoms have 
quantum numbers S = 1/2, L = 0 and 1=1 and 3/2. 

In an optical Feshbach resonance, a laser beam pro¬ 
vides a coupling between the open scattering channel and 
a closed channel containing a bound state dMii ; here 
the open channel is an electronic spin-singlet of ®Li and 
®^Rb, and the bound state is a highly excited LiRb molec¬ 
ular state. When the laser is far detuned from resonance 
with the bound state, the bound state can be adiabati- 
cally eliminated, and we are left with an AC Stark shift 
for the ®Li-®^Rb singlet. The triplet state sees no Stark 
shift. This provides a mechanism for spin exchange. 
While this optically induced spin exchange has not yet 
been experimentally observed, there have been extensive 
studies of both elastic and inelastic scattering proper¬ 
ties near heteronuclear optical Feshbach resonances of 
^Li and ®®Rb Thus the transition frequencies 

for forming ^Li®®Rb molecules are well known. We expect 
that the linewidths, transition matrix elements, and spec¬ 
tral densities for other alkali-metal combinations such as 
®Li®^Rb molecules will be similar. 

Below we provide a mathematical framework to model 
the optical Feshbach resonance and obtain an effective 
interaction between the ®Li and ®^Rb atoms. All the 
physics described in this section is local, and we have 
dropped the index labeling the position of the atoms from 
the second-quantized operators. 

The energy density for the relevant electronic and nu¬ 
clear degrees of freedom in each atom and molecule is of 
the form. 


H = + H^p + Hyaoi + Hpesh- (5) 

.ffmol models the binding energy of the molecule: 

Hraol = ^ ^ I ( 6 ) 

mm' 


where creates a molecule with an electronic spin 

S = 0 and electronic orbital angular momentum J = 1. 
The indices m and m' label the nuclear spins of the ®Li 
and ®^Rb atoms. If the quantization axis of the electronic 
orbital angular momentum is chosen along the direction 
of angular momentum of the laser photon inducing the 
Feshbach resonance, then only one of the molecular states 
in the J = 1 triplet is coupled via the laser to the atomic 
singlet. We denote the binding energy of this molecular 
state by E^. 


The hyperfine Hamiltonians for the atoms are 


ms ,mg 
mi ,m'j 


H^p — hApib ^2 


^(1/2) , 
msm^ 

mimj 

- msmg 

■ 

mimj 


( 7 ) 


ms ,mg 
mi ,m'j 


where h is Planck’s constant, Api = I52MHz and A^b = 
3.4GHz are the hyperfine coupling constants of ®Li and 
®^Rb[2ll, is the vector of spin-S' matrices, and 

Elnsmi ^Insmi Create a ®Li and ®^Rb atom in the 
state \msmi). In terms of the hyperfine eigenstates. 


\ms,mi) = (8) 

F,mp 


where are Clebsch-Gordan coefficients. 

The terms in Elpesh describe the interactions between 
the photoassociation laser and the atoms. We model this 
photoinduced molecular formation by 


Epesh — ^ ^ 


i(fe.r-ajt)-t 

ITT 


- a_ 


V2 


H.c. 


(9) 


where r is the position of the atoms, and Tik and w are the 
momentum and frequency of the laser photon inducing 
molecule formation. The detuning between the atomic 
and molecular states is hco — F{,, and is the transition 
matrix element from the atomic to the molecular state. 

For large detuning, the occupation in the molecular 
state will be small. Therefore we can adiabatically elim¬ 
inate the molecular state and obtain an effective interac¬ 
tion between the ®Li and ®^Rb atoms using second-order 
perturbation theory: 


knt = J2 






\ Eh — hcjj 


V2 




( 10 ) 




Using Eq. (|H]), the operators and can be 

projected into the hyperfine eigenstate basis. Assuming 
that the chemical potential is set such that the F = 3/2 

and F = 2 manifolds are unoccupied, we project iJint 
into the F = 1/2 and F = 1 manifolds. We obtain an 
effective interaction 


Hint ~ ^ ^a^pblibi, ^ ■ 

( 11 ) 

The first term in Eq. m is of the form of Kondo-like 
interactions with 5s = ^ Eb-hw ’ second term a 
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density-density interaction with gn — j eI-Hu: ’ '''^here gs 
and gn were defined in Eq. ©• Generally, in addition 
there would also be intrinsic interactions which modify 
the values of gs and gn in the experiment. To explore 
Kondo physics, gs should be positive. 

If one wanted to exactly produce the Kondo model 
(where gn = 0 ), one could add more photoassociation 
lasers, for example, coupling the electronic spin-triplet 
atomic states. However as we show in Sec. El the pres¬ 
ence of a nonzero does not change the physics. 

A. Experimental and model parameters 

In this section we estimate our model parameters gs 
and gn for a typical experiment performing optical Fesh- 
bach resonance. We also discuss the issue of atom losses 
in optical Feshbach resonances. 

Experiments implementing optical Feshbach reso¬ 
nances typically suffer from high atom loss rates because 
lasers bring the atomic states close to resonance with 
a bound molecular state. The excited molecular states 
have a finite linewidth, and either dissociate into free 
atoms with large kinetic energies or spontaneously de¬ 
cay to ground molecular states. The effect of a finite 
linewidth can be incorporated by making the ac Stark 
shift obtained in Eq. OT complex: 


Eh — hu! + ihj 

The real part of g, Re{g) = ^ ^ 

measure of the interaction strength, and determines 
the magnitude of the model parameters gs and gn- 
The magnitude of the imaginary part of g, Kpx = 

(Eb-hu]y^+{h’y)‘^ ’ inelastic collision rate co¬ 

efficient. 

In experiments in which ^Li and ®®Rb atoms are res¬ 
onantly coupled to a molecular state, the inelastic col¬ 
lision rate co-eflicient typically has a value |Kpa| — 
^ ^ 4 X I0“^^?icm^/s for a moderate laser intensity of 
lOOW/cm^ [ 23 . Typical linewidths are 7 ^ lOMHz. We 
expect that 7 and KpA would have similar values for any 
other alkali-metal combination, and in particular for ®Li 
and ®^Rb as well. We note that is proportional to the 
laser intensity. The inelastic collision rate can be reduced 
by increasing the detuning of the laser. If the laser de¬ 
tuning is 10 times the linewidth {\Eh — fiuj\ = lOTiy), then 
KpA ~ 4 X 10“^^?icm^/s and Re{g) ^ 4 x 10“^^?icm^/s 
for a laser intensity of lOOW/cm^. 

The relevant quantities for estimating the temperature 
scale for observing Kondo physics are and where 
y is the density of fermions. We find in Sec. IIVI that for 
this minimum to occur at a temperature of 0{0.05Te), 
should be 0(0.1). This can be achieved with a 
density of ^ ^ 10 ^^cm“^ and an interaction strength 
l 5 s,n| 2X 10“®?icm^/s, which requires roughly 500 times 


larger intensity than that in [ 2 ^. A judicious choice of 
the resonance may significantly reduce the intensity re¬ 
quired. 

IV. KONDO-ENHANCED SCATTERING 
BETWEEN ®^RB AND ®LI 

Here we calculate the momentum transfer in a colli¬ 
sion between a fermionic cloud and a bosonic cloud. We 
show that spin-exchange collisions lead to a logarithmic 
temperature dependence of the momentum transferred. 
This logarithm is characteristic of the Kondo effect, and 
analogous to the behavior of electrical resistance of mag¬ 
netic alloys. As shown in Fig. [U it leads to a minimum 
in the momentum transferred. The most naive way to 
measure this momentum exchanged would be to launch 
the Bose gas into a stationary Fermi gas and measure the 
final momentum of the Fermi gas. We briefly consider an 
alternative method in Sec. IIV Al 

The duration of interaction between a boson and the 
Fermi gas in the experiment described above is t = L/v 
where L is the size of the Fermi cloud. We calculate the 
momentum transferred from the Bose gas to the Fermi 
gas at time t to zeroth order in \/Mh, first order in v, 
and third order in the interaction parameters gs and gn- 
We perform this calculation for general values of gs and 
gn that are independent of each other. At the end of 
our calculation we specialize to the values of gs and gn 
produced by our proposal in Sec. IHII Since L is a macro¬ 
scopic quantity and we work in the small v limit, we 
make a long time approximation wherever possible. We 
assume that the Bose gas is dilute, and neglect events in¬ 
volving scattering of a fermion with more than one boson. 
Equivalently we calculate the momentum transferred by 
one boson with momentum MhV, and sum over all bosons. 
The Fermi surface will play an important role. 

We consider the collision of the Fermi gas with one 
boson with spin projection m at time 0. The mo¬ 
mentum of the Fermi gas at time t is then = 
( 2 ^ / d^kY,^hknkamit), where 

nham{t) = {bMhV,m{0)cil.^{t)dka{t)b]y^^^^{0)) (13) 

is the occupation of fermions with momentum k and spin 
projection a at time t- In Eq. m the expectation value 
is taken over a thermal ensemble of fermions, with no 
bosons present. The bosonic creation operator preced¬ 
ing the ket state in Eq. m ensures that we calcu¬ 
late the occupation Uk after the collision of one boson 
with the Fermi gas. Since the bosons are spin unpo¬ 
larized, the average momentum imparted by a boson is 

= 3 ( 2^)3 I d^kJ2anib-knkam{t)- Multiplying by Nt, 
the number of bosons, the net momentum of the Fermi 
gas is 

Pit) = / d^k^hknkamit)- (14) 

' / rvm 
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In appendix 13 we describe our diagrammatic pertur¬ 
bation theory approach for calculating rikam (t) ■ We find 
that 

1 f ■ vp{ek) dfk 

q / ^ ^kam\^) —Jk 


y 2 


dck 


fS{S + l)..^ 9^,S{S+1) 

\ A 9s ' 9n 


4(27r)^ 


X / (fp- 


^k 


(15) 


plus terms which scale as or l/Mb- Due to our 

use of point interactions, the interaction parameters ps 
and Qn are renormalized to ps and pn- These renormal¬ 
ized (physical) coupling constants are the ones appear¬ 
ing in Eq. m- This renormalization of the interaction 
strength occurs at all orders of perturbation theory. 

To calculate P{t), we sum the contributions due 
to all momentum states, and include the temperature 
dependence of the fermionic chemical potential, p, = 

ef (^1 - f? (^)^) + ^ 

times, 


^_ 3^(^+l)iV, 

8 



2 

{kpP) hkp 


X 



4 

S{S+1) 




3J 

2ep 


^1.13-b 




+ ilog 


keT 

4eF 




(16) 


where J = ^'^d ^ is the density of fermions. In Eq. 

dni) we have neglected terms which scale as 
or T^. According to our proposal in Sec. IIIIl — 3 

and S' = 1. The result of Eq. (IT6t is plotted in Fig. [T] 
using these parameters and J = O.Icf- For comparison, 
we also plot the momentum transferred to the Fermi gas 
for spin-independent interactions with the same value of 
Pn = —3 X and ps — 0- The logarithmic tempera¬ 
ture dependence of P for spin-dependent interactions is 
characteristic of Kondo physics. Equation (fTHll breaks 
down when log ~ 0(1). Below this temperature, 
the logarithmic increase saturates to a constant. Calcu¬ 
lation of this saturation is the subject of the Kondo prob¬ 
lem and can be addressed with renormalization group or 
Bethe ansatz methods. Equation (TTHl) also breaks down 
when V ~ 

hkp 

The momentum transferred |P| has a minimum at 
a temperature T^m ~ 2 ^ J -P^" 

rameters ^ = —3, 5” = 1, and J = O.Icf, this min- 

93 

imum occurs at a temperature ^ ~ 0(0.05). At 


this temperature and interaction strength, the momen¬ 
tum imparted by one boson to the Fermi gas is ~ 

jhkp [kpL). For a 20^m long Fermi cloud at a 

density of 10^^cm“^, the momentum imparted per boson 
is nearly 1.2fikF- We estimated in Sec. [TTIlthat achiev¬ 
ing J = 0.1 ef would require high intensity lasers and 
tight trapping of the fermions. The observation of this 
minimum will be a direct experimental confirmation of 
Kondo physics. 


A. Alternative methods to measure enhanced 
Kondo scattering 

Here we briefly explain an alternative method to mea¬ 
sure the enhanced Kondo scattering between a Fermi 
cloud and a Bose cloud. We consider inducing dipole 
oscillations of a Bose cloud and a Fermi cloud in a har¬ 
monic trap of frequency w. The clouds will collide ev¬ 
ery half-cycle and exchange momentum P. As a result 
the amplitude of oscillations of the Fermi cloud will re¬ 
duce each half cycle. Conservation of momentum implies 
that the maximum fermion displacement X will reduce 

each half cycle by 5X ~ ^ where Na is the number 

of fermions; the Bose clou3’s amplitude will not change 
very much because of the bosons’ heavy mass. The Bose- 
Fermi interaction interval is longer for a smaller relative 
momentum, and vice versa. Thus the momentum ex¬ 
changed I PI is independent of the relative velocities of the 
cloud, leading to a linear decay of the amplitude rather 

than exponential; ^ ^^® Bose-Fermi in¬ 

teractions are Kondo-like, the damping rate of amplitude 
of oscillations will have a minimum at the same tempera¬ 
ture as I PI does, T^in I '^%‘i '■ a typical 

V ^''"§fs(s+i) 

amplitude of oscillation X ~ lOO^m in a trap of fre¬ 
quency ct) = 27r X 10 Hz, and if ^ = 200, the amplitude 
will decay to zero in about 12 oscillations at T = Tmin- 
The observation of a minimum in the damping rate will 
also be an experimental confirmation of Kondo physics. 


V. SUMMARY 

We considered scattering between a spin-1/2 Fermi 
gas and a dilute spin-unpolarized Bose gas. As an 
example we considered ®Li and ®^Rb as our itinerant 
fermions and bosonic magnetic impurities. We proposed 
using an optical Feshbach resonance to produce rotation- 
ally symmetric interactions between the ®Li and ®^Rb 
atoms, which included both spin-dependent Kondo-like 
and spin-independent density-density interactions. We 
argued that these interactions would give rise to en¬ 
hanced Fermi-Bose scattering. We perturbatively cal¬ 
culated the temperature dependence of the momentum 
transferred to the Fermi gas in a scattering experiment. 
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up to third order in the Bose-Fermi interaction strength. 
We showed that the temperature dependence of the mo¬ 
mentum transferred has a minimum at a characteris¬ 
tic temperature and is logarithmic at low temperatures, 
characteristic of the Kondo effect and analogous to the 
behavior of electrical resistance in magnetic alloys. 

Our proposal to implement spin-dependent interac¬ 
tions requires overcoming significant experimental chal¬ 
lenges such as using high intensity lasers to achieve large 
interaction strengths. However, overcoming these chal¬ 
lenges enable the possibility of exploring exotic phe¬ 
nomena due to Kondo physics. The ground state of a 
Bose-Fermi mixture with Kondo-type spin-dependent in¬ 
teractions should display interesting correlations, with 
each boson surrounded by a screening cloud of fermions 
with opposite spin Q. These clouds may be observable 
through various imaging techniques [22l - l24| . Similar ex¬ 
periments with bosons confined to a lattice would probe 
an analog of the Kondo lattice problem. 

One can explore other techniques to experimentally 
produce Kondo-type interactions. For example, optically 
coupling the electronic triplet states of ®Li-®^Rb with ex¬ 
cited molecular states will lead to a rotationally asym¬ 
metric interaction which also displays Kondo physics. 
Alternatively, one can realize the Anderson model and 
Kondo-like situations by trapping impurities in deep po¬ 
tentials [25l - [^ . 
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Appendix A: CALCULATION OF THE 
MOMENTUM TRANSFERRED 


Here we calculate nkam{t) in Eq. (IT^ and P{t) in 
Eq. (fT4l) . The standard way to calculate quantities like 
nkam{t) is using the S'-matrix [^ : 


v,m (0)4a(Oa/ca(O^L,«,m(0))o, 
g — g-i f dTHi^ir) ^ 


(Al) 


where T orders the operators along a path shown in Fig. 

which starts at time 0, passes through time t, and re¬ 
turns to time 0. All our integrals over time follow this 
path. The notation ()o implies that all operators inside 
()o evolve according to 


afca(t) = 

bk^it) = 


In Eq. (IA3I) . 0(ti — ^ 2 ) = 1 if ti is after t 2 along the path 
in Fig. [21 and 0 otherwise. 

We perturbatively expand nkam{t) in the vertex de¬ 
picted in Fig. ISDc), whose value is 



+ ^3 - ^2 - ^4) 

(A4) 


The vertex denotes a scattering event between a 
fermion and a boson. The time at which this scattering 
event occurs is integrated over the path in Fig. [2] All 
momenta and spin projections are summed/integrated 
over, with the constraint that momenta and spin are con¬ 
served at each vertex. The diagrams which contribute to 
Eq. (ED have four external propagators. There is an 
incoming and outgoing fermion propagator evaluated at 
time t, and carrying momentum hk and spin projection 
a. There is also an incoming and outgoing boson propa¬ 
gator evaluated at time 0, and carrying momentum MbV 
and spin projection m. All lines and vertices in a Eeyn- 
man diagram can be labeled using the rules described 
above. Therefore we omit labels. Einally, each diagram 
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(b) 

FIG. 3: Diagrammatic representation of vertex and propa¬ 
gators. (a) Solid line denotes a fermion propagator which 
propagates a fermion with momentum k and spin projection 
a from time t 2 to ti. (b) Dashed line denotes a boson propa¬ 
gator which propagates a boson with momentum k and spin 
projection ^ from time t 2 to ti. (c) A vertex denotes the ma¬ 
trix element for a Bose-Fermi scattering event. Mathematical 
expressions are given in Eas. (IA3|l and (IX4l) . 



FIG. 5: First-order diagrams in the expansion of nkam{t). 



FIG. 6: Two of the diagrams that are zero at second order. 


FIG. 4: Zeroth-order diagram in the expansion for nkam(t)- 


carries a multiplicity, which is the number of times it ap¬ 
pears in the expansion of Eq. m in powers of gs and 
9n- 


2. Calculation of nkam(t) 

Terms of in the perturbative expansion of 

nkamit) contain 2n -|- 2 pairs of operators leading to 
(n-|-1)!^ contractions. The resulting number of diagrams 
increases exponentially with n. We explicitly consider 
each order and evaluate the nonzero diagrams. 


a. Zeroth order 

The expression for the zeroth-order term in the expan¬ 
sion of nkam{t) is 

(A5) 

Using Wick’s theorem, 

= (^M,v,m(0)6L,«.m(0))o(aLWafc«(t))o 


= fk. 


(A6) 


The corresponding Feynman diagram is shown in Fig. 
m Since the bosons and fermions do not interact at this 
order, rSklm doss not contribute to any momentum trans¬ 
fer. 


b. First order 

The first-order term in the expansion for Ukamit) is 


kam 


[t) = -i j cLti 


(A7) 

By Wick-contracting the above expression, we find that 
is the sum of the four diagrams shown in Fig. [SJ 
all of which evaluate to zero. For example, 




= y dT = 0. (A8) 


Due to the same reason. Figs. [5jb), [5jc) and [5jd) are 
also zero. Therefore, 



1 


= 0 - 


(A9) 


Moreover, the same reasoning implies that all higher- 
order diagrams in which a fermion or boson loop begins 
and ends at the same vertex are also zero. 


c. Second order 


The second-order term. 


/ dTidT2(T7fint(Ti)7fint(T2)&Mt,;,m(0) 




(AlO) 











(a) 



(b) 



FIG. 7: Nonzero diagrams at O(g^) in the expansion for 
(t) ■ 


can be contracted into Wick pairs in 36 ways, which give 
rise to 20 different diagrams. Most of these diagrams 
are zero because of reasons explained in Sec. lA 2 bl In 
addition, the diagrams shown in Fig. [5] also evaluate to 
zero. For example, since we work in the dilute boson 
limit, there can only be one boson line in any time slice, 
implying that Fig. IHta) is zero. The only two nonzero 
diagrams are shown in Fig. [3 
Using our Feynman rules, 




V(2TTf 


d^P (1 - fk)fi 


sin^ 5et/Ti 


<5e2 


2S{S+1) ^ 2 

X I g. . ' + 2 g^ 


and 


1 ^ 
^ 2-^ 


V{2t:Y 


J (Pp /fc(I - /p)- 




(All) 


sin^ Se'tjh 
(A12) 


where 5e = 5 (ee - e,, - 

S/ = l(e,-e, + lAfe.^- <«‘-« 

glecting terms of order l/M^, 6e = Se' = 
1 {ek-m^v/h - ep-mav/h) ■ The resulting second-order 
contribution is 


:E 


kotm 


it) = - 


„ 2 S(S+ 1 ) 
ys 2 


2 g^ 


U(27r)3 


J d^p{fk-fp) 


sin^ t6e/Ti 
(5e2 


(A13) 

Since the bosons are much heavier than the fermions, 
they have nearly the same velocity v before and after 
scattering. Therefore, it is easier to work in the bosons’ 
rest frame. For small v, 


it)=- 

o / ^ k-\—' 


9. 


2S(S-H) 


+ 2 g^ 


d^P 


U(27r)3 

X \ fk - fp + ‘hk ■ - Tip - V 


dfk -.dfp 


dek 


dir. 




((Cfc - ep)/ 2)2 


(A14) 


where 0{v'^ ,1/Mb) refers to terms which scale as or 
l/Mb- The first two terms in Eq. (IAI4|) have negligible 
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FIG. 8: Non-zero diagrams at 0{g^) in the expansion for 
‘kl'kam (t). 


contribution near Cfc = Cp. At long times, any signifi¬ 
cant contribution comes from the tail of ) 

where sin^(t(efc — ep)/2fi) can be approximated by its 
average, 1/2. Hence their contribution saturates to a 
constant at long times. For the last two terms in Eq. 
(IA14I) . which are significant near = Cp, we approxi¬ 
mate Hence at long times, 


1 ( 2 ) , , _ + 2 g^)t 


{t) 

/ J k-\—' 


U( 27 r) 


(A15) 


.9/fc 




= -4 


dek 


2 S{S+1) 2 

ys 4 ' yn 


de. 


(Tp k ■ - p- V— 6{ek -ep) + 0[ f, v^, — 


Mb 


1/2 


tk ■ v^^p{ek) +0 [t^, v^, 


de 


Mb 


where p{ek) is the three-dimensional density of states for 
a single spin projection. In the laboratory frame. 




-4 


„ 2 S(S+l) 
ys 4 


1/2 


-tk 


-.dfk , , 
v^p[ek) 
dek 


+ 0 



(A16) 


d. Third order 
The third-order term 

yimi*) j dTidT2dT3{THint{Ti)x (Al7) 

-^int (^2 )-^int (^3 (0) “la {k)dka (0))o 

can be contracted into Wick pairs in 576 ways. However 
due to reasons explained in Secs. lA 2 bl and lA 2 cl all dia¬ 
grams except the ones shown in Fig. [5] are zero. After a 



































treatment similar to the one at second order, we calculate 
the third-order contribution to be 


9 




-Sg^^gnSiS + 1) - 4gl) + O 

(A18) 


The right hand side of Eq. (IA18I) consists of an ultra¬ 
violet divergent term arising from J , and a fi¬ 
nite term f which will ultimately give rise to 

a logarithmic temperature dependence. The ultraviolet 
divergence is an artefact of choosing a contact potential 
between the fermions and bosons which is nonzero only 
when they are at the same location in space. In reality, 
the interaction between the fermions and bosons has a 
finite range, which removes the ultraviolet divergence by 
introducing an upper cutoff on the limits on the integral 
over momenta. The exact details are unimportant if we 
express our results in terms of physical quantities. To 
this effect, we define effective coupling constants gs and 
gn where 


~2 2 
9 s = 9s 


~2 2 
9n=9n 


1 -^ 


1 -h 


9s + 6 g„ 
2(27r)3 

9n / 


d^p- 


1 


(27r)3 


dV- 


1 


Cfc ~ Cp 


(A19) 


The result for nkam{t) has no ultraviolet divergences 
when expressed in terms of gs and gn- 
The resulting nkam{t) at long times is 


1 

3 


m 


fk 


4tk-vp{ek) dfk 
E 2 dek 


f SjS + l) ^.^ g^,S{S + l) 

V 4 4(27r)3 


J 

(A20) 


3. Final momentum of the Fermi gas 

The total momentum P of the Fermi gas [defined in 
Eq. (IT4l) ] will be along the direction of v. Its magnitude 
is 



SthNi, 
vV (27r)3 



dfk , - 


5(5- 


—fs + gl 


-glsjspi) 

4(27r)3 



(A21) 


After integrating out the angular co-ordinates of k and p 
and performing a change of variables, 


|P| =- 


IdrUaLNi, 

3hV^ 




de 


' 5 ( 5 + 1 ) .2 -2 
X I —-gl + gl 


glsjs + i) 

AV 


de- 


P(ep)/(ep) 


(A22) 


We evaluate the second-order terms using a Sommer- 
held expansion. 


17^2 I 


4?ieF 


+ I)(l + a2) 



keT 

cf 



(A23) 


where J = gsy and a 


Qn _ 2 _ 

9s y/S{S+l)' 


The third-order terms are 




45(5 + l)maLNb f g. 


Sh 


df{e) 


V “‘‘as 


X J 


de, 


95(5 + l)m,LW ^3 r ^2 df{e) 


16Tie 


9/2 


de 


de- 


pV'^p 


P(ep)/(g 

^ f _ f 

C Cp 

/(ep) 

e- en' 


(A24) 
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We simplify the above expression by performing integration by parts, 

3, 9S(S+l)m.LJVs ,3 j3/(e) f , S f ^ 

L 

F 


JQ 


/O 


Cp 


9S'(S' + l)maLNb 


She 


9/2 




/O 


5e 


/o 


cp 


2 


log 




^ log 


•\/e + 


/3(\/e + 


We split Eq. (IA25I) into two terms. We evaluate one of these terms numerically, 


o 


ksT Y 

CF j 


We use a Sommerfield expansion for the remaining term. The result is 
9S{S + l)maTiVb 


1^31=^-- 


She^p 




(A25) 


(A26) 


48/ y ep 


O 


ksT 

ef 


The final momentum of the Fermi gas is 

P =Po^ f 1 + ^ f ^ V - f 1.13 T f 2.6 - ^ 

I 6 \ ep J 2(l + Q!2)e_F \ \ 48 


(A27) 


{ksTV 1 , kpT ^ 5 f 2 fkBTY 


(A28) 


where Pq = ^ and as before, we neglect terms of O (t°, v^, • 


[1] J. Kondo, Prog. Theor. Phys. 32, 37 (1964). 

[2] A. C. Hewson, The Kondo problem to heavy fermions, 
Vol. 2 (Cambridge university press, 1997). 

[3] I. Affleck, arXiv:0911.2209 (2009). 

[4] I. Affleck, in Strongly Correlated Fermions and Bosons 
in Low-Dimensional Disordered Systems (Springer, New 
York, 2002), pp. 1-12. 

[5] J. B. Boyce and C. P. Slichter, Phys. Rev. B 13, 379 
(1976). 

[6] J. B. Boyce and C. P. Slichter, Phys. Rev. Lett. 32, 61 
(1974). 

[7] H. Tsunetsugu, M. Sigrist, and K. Ueda, Rev. Mod. Phys. 
69, 809 (1997). 

[8] Z. Fisk, J. L. Sarrao, J. L. Smith, and J. D. Thompson, 
Proc. Natl. Acad. Sci. 92, 6663 (1995). 

[9] Z. Fisk, H. Ott, and G. Aeppli, Jap. J. Appl. Phys. 26, 
1882 (1987). 

[10] R. A. Williams, L. J. LeBlanc, K. Jimenez-Garcia, M. C. 
Beeler, A. R. Perry, W. D. Phillips, and I. B. Spielman, 
Science 335, 314 (2012). 

[11] N. R. Thomas, N. Kjmrgaard, P. S. Julienne, and A. C. 
Wilson, Phys. Rev. Lett. 93, 173201 (2004). 

[12] S. D. Gensemer, R. B. Martin-Wells, A. W. Bennett, and 
K. Gibble, Phys. Rev. Lett. 109, 263201 (2012). 

[13] N. Robins, G. Figl, J. Close, et ah. Optics express 16, 
13893 (2008). 

[14] H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. 
Rev. Lett. 91, 250402 (2003). 

[15] S. Dutta, J. Lorenz, A. Altaf, D. S. Elliott, and Y. P. 


Chen, Phys. Rev. A 89, 020702 (2014). 

[16] P. O. Fedichev, Y. Kagan, G. V. Shlyapnikov, and 
J. T. M. Walraven, Phys. Rev. Lett. 77, 2913 (1996). 

[17] K. M. Jones, E. Tiesinga, P. D. Lett, and P. S. Julienne, 
Rev. Mod. Phys. 78, 483 (2006). 

[18] C. Ghin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 
Mod. Phys. 82, 1225 (2010). 

[19] M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, 
G. Ruff, R. Grimm, and J. H. Denschlag, Phys. Rev. 
Lett. 93, 123001 (2004). 

[20] S. Dutta, D. S. Elliott, and Y. P. Chen, EPL (Euro¬ 
physics Letters) 104, 63001 (2013). 

[21] E. Arimondo, M. Inguscio, and P. Violino, Rev. Mod. 
Phys. 49, 31 (1977). 

[22] L. W. Cheuk, M. A. Nichols, M. Okan, T. Gersdorf, V. V. 
Ramasesh, W. S. Bakr, T. Lompe, and M. W. Zwierlein, 
Phys. Rev. Lett. 114, 193001 (2015). 

[23] M. F. Parsons, F. Huber, A. Mazurenko, C. S. 
Chiu, W. Setiawan, K. Wooley-Brown, S. Blatt, and 
M. Greiner, Phys. Rev. Lett. 114, 213002 (2015). 

[24] M. Miranda, R. Inoue, Y. Okuyama, A. Nakamoto, and 
M. Kozuma, Phys. Rev. A 91, 063414 (2015). 

[25] M. Foss-Feig, M. Hermele, and A. M. Rey, Phys. Rev. A 
81, 051603 (2010). 

[26] J. Bauer, G. Salomon, and E. Dernier, Phys. Rev. Lett. 
Ill, 215304 (2013). 

[27] Y. Nishida, Phys. Rev. Lett. Ill, 135301 (2013). 

[28] M. Nakagawa and N. Kawakami, Phys. Rev. Lett. 115, 
165303 (2015). 
































11 


[29] 


G. D. Mahan, Many-particle Physics (Springer, New 
York, 2000). 



